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Abstract 

We investigate the e + p — > fi + X event with high transverse momenta observed 
at HERA (HI) and show that this event could be interpreted as a signature 
of the single production of the scalar top quark in a supersymmetric model 
with imparity breaking interactions. The event topology of the HI event is 
rather characteristic and in fact it can be simulated by our specific scenario if 
we reasonably choose our model parameters to be (i) tutt- ~ ITeV [0.8TeV] 

for A' 131 = 0.1 [0.05] and (ii) ~ 150GeV, lOOGeV ~ m~ ~ 200GeV and 

A' 131 ~ 0.05. 
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Recently the HI group reported that an event with high transverse momenta had been 
observed [0 at HERA. The total data sample analyzed at the HI group corresponds to 
an accumulated luminosity of 3.2pb _1 in positron (27.5GeV) - proton (820GeV) collisions 
and of O.Spb^ 1 in electron - proton collisions. The event is characterized by (i) a large 
transverse momentum of the single muon fi + , Pt{p) = 23.4 db 2.4^t 5 GeV, (ii) a small 8 
(= ~ cos0)), 5 = 19.2 ± 1.6l2jGeV, where E and 6 denote the energy and angle 

of any detected particle, (hi) a large transverse momentum of the total hadronic system, 
-Pr(hadron) = 42.1 ± 4.2GeV and (iv) a large missing transverse momentum Pr(miss) = 
18.7 ±4.8lf GeV. 

Some possible interpretations of the event have been given in ref.Q; 

(A) production of high Pt jets 

(B) W production and its leptonic decay 

(C) A flavour changing neutral currents (FCNC) or leptoquark production 

However, it does not seem to give reasonable explanation of the event characterized by the 
topology mentioned above Jl|. As for the scenario A, the probability for the event being 
due to the production of two high Pt jets, where one jet shows the signature of a muon, 
is smaller than 10 -3 . As for the scenario B, the W production and its leptonic decay 
would give a rather small transverse momentum of the total hadronic system. In fact, the 
Monte Carlo calculation shows that for Pr(hadron) > 40GeV the cross section is reduced 
to 7fb. That is, with one event seen in 4pb~ 1 we are left with a 3% probability for this 
interpretation of the event. As for the scenario C, expected events would have balanced 
Pt and a value of 5 = 2E e . This is due to the fact that events originated from the FCNC 
or leptoquark production would show topologies like neutral current deep inelastic events, 
but with the final state positron replaced by a muon. The fact that less than 1% of neutral 
current events show a value of 5 < 20 GeV disfavours this interpretation. 

In this letter, we propose a possible explanation of the single muon event in the frame- 
work of the minimal supersymmetric (SUSY) standard model (MSSM). We will show that 
the high Pt muon could appear from the single scalar top quark (stop) production through 
an i?-parity breaking coupling at HERA. 

In the previous work Q, we have already shown that one of the signals of the single 
stop production to be detected at HERA is characterized by the high Pt spectrum of 
muons. First, for the sake of convenience we will briefly summarize the basic idea by 
referring to our previous work [^. The discussion is based on the MSSM with an i?-parity 
breaking (RB) interaction 

L = A' 131 cos e t (hdP L e + t\eP R d), (1) 

where A' 131 and 9t, respectively, denote the coupling strength and the mixing angle of 
the stops ||, |J. Here Pl,r read left and right handed chiral projection operators. The 
interaction Lagrangian (|l|) has been originated from the general RB superpotential M; 

W H = XijkLiLjE-k + X' ijk LiQjD c k + X'^U^D^D^, (2) 
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where i,j,k are generation indices. The first two terms violate the lepton number L 
and the last term violates the baryon number B. If we want to explain such unresolved 
problems as (i) the cosmic baryon number violation, (ii) the origin of the masses and the 
magnetic moments of neutrinos and (iii) some interesting rare processes in terms of the L 
and/or B violation, the i?-parity breaking terms must be incorporated in the MSSM. 

The coupling Eq. (jlj) will be most suitable for the ep collider experiments at HERA 
because the stop will be produced in the s-channel in e-q sub-processes [||, ||] 

ep -> t x X. (3) 

Note that the stop cannot couple to any neutrinos via i?-breaking interactions. This is a 
unique property of the stop which could be useful for us to distinguish the stop from some 
leptoquarks. Production processes of the first and second generation squarks have been 
discussed in ref.j/J. 

Babu and Mohapatra H have recently shown that the severe constraint on a product 
A' 113 A / 131 ~ 3 x 1CP 8 comes from experimental data of the neutrinoless double (5 decays. 
Here, we will assume A' 131 to be only non-zero coupling parameter in what follows. The 
upper bound on the strength of the coupling has been investigated through the low-energy 
experiments [[| and the neutrino physics H. The most stringent bound A' 131 ~ 0.25 comes 
from the atomic parity violation experiment ||. 

Next we examine the decay modes of the stop. In the MSSM, the stop lighter than 
the other squarks and gluino can decay into the various final states : 
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where (k = 1 ~ 4), Wi(i = 1,2), v and £, respectively, denote the neutralino, the 
chargino, the sneutrino and the charged slepton. (a) ~ (g) are the i?-parity conserving 
decay modes, while (h) is only realized through the RB couplings (|l|). 

If we consider the stop with mass small enough in the case of the R conserving coupling, 
the first five decay modes (a) to (e) are kinematically forbidden due to the observed top 
mass mt— 175 GeV [l(| as well as the model independent lower mass bounds for sparticles 
; m^y ~ 45 GeV, mj ~ 45 GeV and m~ ~ 40 GeV. So (f) and (g) survive. Hikasa and 

Kobayashi Q have shown that the one-loop mode (g) t\ — > cZ\ dominates over the four- 
body mode (f) ti — > bff'Z\. So we can safely conclude that such a light stop will decay 
into the charm quark jet plus the missing momentum taken away by the neutralino with 
almost 100% branching ratio. On the other hand, if we consider the RB coupling A' 131 > 
0.01, which roughly corresponds to the coupling strength to be detectable at HERA, the 
decay modes (c) to (g) are negligible due to their large power of a arising from multiparticle 



3 



final state or one loop contribution. Then only two body decay modes (a), (b) and (h) are 
left for our purpose. 

We have found Q that if the stop is heavy enough, i.e., > m b + m^> and the RB 

coupling is comparable with the gauge or Yukawa coupling X' ^/Att ~ a, at there is a 
wide range of parameters where BR(t\ — > bW^) dominance over BR(t\ — > e<i) is assumed. 

In this case, we should take into account of the process 

ep -» 6^ fc X, (4) 

where the virtual contributions of the sneutrino with the same RB coupling constants A' 131 
have, of course, been considered. The differential cross section is given by 

da , ~ _ oAl r .„ (« " "»&)(« " m L ) 
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with s = xs,i = —Q 2 and 

G L ee - m ^ 2C ° 5 ^ , (6) 
V2mvK cos (3 

G R = Vfci cos 6> t + -= — . (7) 

\/2mw snip 

Here Vki and stand for the chargino mixing angles |ll| . The mixing angles as well 



as masses of the neutralinos m-x and the charginos m~, are determined from the basic 
parameters in the MSSM (//, tan/?, M%). We can see that the e + beam is more efficient 
than the e~ one to distinguish the stop signal from the SM background. This can be 
understood from the fact that the e + collides with valence d-quark in the proton, while 
the e~ does only with sea d-quarks. It is expected that the detectable cross sections 
a ~ 0.1 pb for heavy stop with mass m~ ~ 250 GeV for e + beams. As far as e~ beams 

are concerned e~p — ► bW^X would be detectable for m~ ~ 170 GeV. In our model the 

LSP, the lightest neutralino Z\ possibly decays into i?-even particles via only non-zero RB 
coupling A' 131 . A typical decay chain will be 

ep -» bWiX -» {UvZ x )X -» b{lv(bdv))X. (8) 



Thus, a possible typical signature of the stop production ep — > bW\X would be b- 
jet+lepton+fV in the case of no LSP decay or 26-jets+jet+lepton+ in the case of the 
LSP decay via RB coupling. One of the signals to be detected at HERA is characterized 
by the high Pt spectrum of muons. The lower Pt cut certainly makes the event distinctive 
from its background. The cross section a(e + p — * t\X — > bW\X) varies from 1 to lOpb 
depending on mass of the stop in the range of 100 ~ 150GeV. 
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Now it is the position to present our calculation for some kinematical distributions in 
the process (||), which will also be compared to the experimental distributions of the HI 
event ]|]. 

First, the Pt{ij) distribution of the expected number of events is shown in Fig.l. In 
the calculation, we take a typical set of model parameters, (fi, M2, tan/3, m t , 9t, A' 131 ) 
= (-300GeV, 50GeV, 2, 175GeV, l.Orad, 0.1) and the integrated luminosity S^pb" 1 . In 
this case we get the lighter chargino mass = 63GeV and the lightest neutralino mass 

to- = 28GeV. For simplicity, the branching ratio BR{W\ —* v[iZ\) is assumed to be 
5 EMl ■ ^ ne dependence on the branching ratio of the chargino will be discussed later. We 
find in Fig.l that rather heavy stop, mr = 100 ~ 120, could give a high Pt{^) event at 
the present integrated luminosity. 

We show the 5 distribution together with the experimental data in Fig. 2. The observed 
value of 5 is significantly smaller than the allowed maximum value 2E e . Such 5 could be 
obtained only when the LSP does not decay via the RB couplings in the detector. This 
constraint leads to following important consequences. We can quantitatively present the 
requirement as cir^ ~ lm, which corresponds to T~ ~ 10 eV. By calculating the 3- 

body LSP decay width, we get severe constraints, ~ ITeV [0.8TeV] for A' 131 = 0.1 

[0.05]. Another one is that the total transverse hadronic momentum should be supplied 
by a b quark at the first vertex in the process (g). In other words, we can take Prib) 
~ P^(hadron) in our model calculation. Note, moreover, that large masses of sfermions 
justify our assumption on the chargino decay, BR(W\ — ► vpiZ\) ~ g. When sfermions are 
sufficiently heavy the dominant contribution to the decay matrix elements comes from the 
H^-boson exchange diagram, which has apparently universal fermion couplings. 

Shown in Fig. 3 is the Prip) (— Pr(hadron)) distribution. From the figure, we find that 
heavier stop is favourable to simulate the large transverse hadronic momentum ~ 40GeV. 
Note that the maximum value of Pr(b) is determined from the simple kinematics of the 
two-body decay, 

Consequently, we can restrict possible regions in (m^ , to^) parameter space in terms of 
the experimental constraint, Pr(b) ~ Pj. a:p (hadron) ~ 40GeV. 

We present such a experimentally favourable region (shaded area) in Fig. 4. Horizontal 
(dotted) lines correspond to one event with Pt{^) > 20GeV to be expected. It is seen 
from the figure that constraints have been set on the masses and , if we seriously 
take into account the present data. For example, the chargino mass is smaller than about 
lOOGeV for A' 131 ~ 0.1 under the condition of Pr{b) ~ 40GeV. Moreover, A' 131 should be 
larger than about 0.05 and the stop mass should be larger than about lOOGeV because of 
the LEP bound on to^ > 45GeV. As has already been mentioned, low energy experiments 

set the constraint A' 131 ~ 0.25. 

We have investigated a possible scenario to explain the HI single muon event by the 
single production of the stop with an i?-parity breaking interaction in the framework of 
the MSSM. The event topology of the HI event is rather characteristic and in fact they 
can be simulated by our specific scenario if we restrict arbitrary model parameters to some 
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reasonable ranges. In order to work our scenario, we must have (i) fnji~ ~ ITeV [0.8TeV] 
for A' 131 = 0.1 [0.05], (ii) ~ 150GeV, lOOGeV ~ m~ ~ 200GeV and A' 131 ~ 0.05. 

Our scenario would be confirmed or rejected at LEP2 or next linear colliders through 
the search for e + e~ — > W^W^. Certainly, the discovery of the chargino could reveal us 
that the nature is supersymmetric but could give no information as to whether or not the 
nature does not respect the i?-parity. To confirm it, we should seek for the stop with the 
i?-parity breaking interaction, e.g., through a search for e + e~ — > edt\. Needless to say, for 
our purpose it would be highly desirable to carry out the high luminosity run at HERA. 
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Figure Captions 



Figure 1: Pt(i^) distribution of the expected number of events together with the experi- 
mental data. We take (//, M 2 , tan/3, m t , 6 t , A' 131 ) = (-300GeV, 50GeV, 2, 175GeV, 
l.Orad, 0.1) and the integrated luminosity as S^pb^ 1 . Solid line and dotted line 
respectively correspond to m~ =120 GeV and m~ =100 GeV. 

Figure 2: 5 distribution of the expected number of events together with the experimental 
data. We take mj =100 GeV and the same parameters in Fig.l. Solid line and 
dotted line respectively correspond to no LSP decay and the LSP decay within the 
detector. 

Figure 3: Pr(b) (— ^(hadron)) distribution of the expected number of events together 
with the experimental data. Parameters are the same in Fig.l. 

Figure 4: Favourable region (shaded area) in (m^, rrij^) parameter space by the experi- 
mental constraint, Pr{b) ^ P^^hadron) = 42.1±4.2GeV. Horizontal (dotted) lines 
correspond to one event with Pt(/j) > 20GeV to be expected. 
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